Electron correlation and Fermi surface topology of Na2;Co02 
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The electronic structure of Naa;Co02 revealed by recent photoemission experiments shows im- 
portant deviations from band theory predictions. The six small Fermi surface pockets predicted by 
LDA calculations have not been observed as the associated e'g band fails to cross the Fermi level for 
a wide range of sodium doping concentration x. In addition, significant bandwidth renormalizations 
of the t2g complex have been observed. We show that these discrepancies are due to strong electronic 
correlations by studying the multi-orbital Hubbard model in the Hartree-Fock and strong-coupling 
Gutzwiller approximation. The quasiparticle dispersion and the Fermi surface topology obtained in 
the presence of strong local Coulomb repulsion are in good agreement with experiments. 

PACS numbers: 71.27.-|-a, 71.18.-fy, 74.25.Jb, 74.70.-b 
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The cobaltates (Na2:Co02) are doped 3d transition 
metal oxides in which the Co atoms form a layered 
hexagonal lattice structure. In contrast to the high- 
Tc cuprates, where the Cu^+ has a 3d^ configuration 
and occupies the highest single eg [(1^2 _y2) orbital near 
the Fermi level, the cobaltates are multi- orbital systems 
where the 00^+ is in the 3d^ configuration, occupy- 
ing the three lower t2g orbitals, similar to the ruthen- 
ates (Sr2Ru04). The unexpected discovery of a su- 
perconducting phase of yet unknown origin in hydrated 
Naa;Co02 around x ~ 0.3 has generated renewed inter- 
ests in this material. However, such basic issues as the 
low energy electronic structure and Fermi surface topol- 
ogy in the cobaltates have not been well understood. 
Band structure (LDA) calculations find that the trig- 
onal symmetry of the Co site in the triangular lattice 
splits the three t2g complex into an aig and two degen- 
erate e'g states at the zone center (F point). The LDA 
predicts a large Fermi surface (FS) associated with the 
aig band enclosing the F point and six small FS pockets 
of mostly e'g character near the K points 0, Q . 

However, recent angle-resolved photoelectron spec- 
troscopy (ARPES) measurements on the cobaltates re- 
vealed only a single hole-like FS centered around the F 
point for a wide range of Na concentration x 0, 0, S ■ 
The area enclosed by the FS exhausts the Luttinger vol- 
ume, which is consistent with the observation that the 
dispersion of the e'g band associated with the LDA FS 
pockets lies below and never crosses the Fermi level 
The absence of the FS pockets is unexpected and puts 
serious constraints on several proposed theories of non- 
phonon mediated superconductivity as well as magnetic 
properties based on the nesting conditions of the FS 
pockets 00, El- Furthermore, the measured quasipar- 
ticle bandwidths are significantly narrower than the LDA 
predictions These fundamental discrepancies be- 

tween ARPES and LDA suggest that the effects of strong 
electronic correlations are important in the cobaltates. 
The effects of local Coulomb repulsion U has been con- 



sidered in the LSDA-I-U approach, which indeed finds the 
absence of the small FS pockets However, the lat- 
ter is tied to the fully polarized ferromagnetic state in 
the LSDA-t-U theory which gives spin-split bands and a 
spin polarized FS with an area twice as large. This is 
inconsistent with ARPES and likely an artifact of the 
LSDA-I-U approximation. A recent calculation based 
on the multi-orbital Hubbard model and the dynamical 
mean-field theory finds that the FS pockets become even 
larger in size than the LDA predications [l^ . 

The focus of the present work is to explain how strong 
correlations drive orbital polarization and the band nar- 
rowing observed in ARPES. We adopt a multi-orbital 
Hubbard model description where the nonintcracting 
part is determined by fitting the LDA band structure. 
The interacting part contains both the intra (U) and 
the inter-orbital ([/') local Coulomb repulsion as well as 
the Hund's rule coupling Jh- First, a basis independent 
Hartree-Fock (HE) calculation is performed which is in 
essence a LDA-I-U calculation in the paramagnetic phase. 
We find that for U' much less than U, multi-orbital oc- 
cupation is favored in order to reduce the cost of double 
occupation. As a result, the HE self-energy renormalizes 
the atomic level spacing in such a way that the size of 
the FS pockets associated with the e'g band grows. This 
trend is however reversed when U' grows and becomes 
comparable to U. In the HE theory, the size of the e'g FS 
pockets begins to shrink for U'/U > 3/5. To correctly 
capture the physics of strong correlation for large U and 
U' , we generalize the Gutzwiller approximation to the 
case of multi-orbitals. We find that in the strong-coupling 
regime, orbital polarization is tied to Pauli-blocking, i.e. 
the orbital occupation dependence of the Gutzwiller band 
renormalization factors. We obtain both band narrowing 
and the disappearance of the FS pockets in good agree- 
ment with the ARPES experiments. 

We start with the multi-orbital tight-binding model on 
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a two-dimensional triangular lattice, 

ij,cr af3 i.iy 

(1) 

where the operator rfj^j^ creates an electron in the a or- 
bital with spin cr on the Co site and iij,Q/3 is the hopping 
integral between the a orbital on site i and the /3 or- 
bital on site j. The relevant valence bands near the FS 
consist of the Co t2g = {dxy, dyz, d^^} orbitals and have 
an electron occupancy of 5 -I- a;. The A in Eq. de- 
scribes the crystal field due to trigonal distortion that 
splits the t2g complex into a lower aig singlet and a 
higher e'g doublet, where aig = {dxy + dyz + dzx)/V5, and 

= {{dzx-dyz)/V^, {2dxy-dyz-dzx)/VG}- For conve- 
nience, we will work in the hole-picture via a particle-hole 
transformation d —> d'' , in which the band filling of holes 
is 1 — X. The structure of the tight-binding Hamiltonian 
in k-space on the triangular lattice is 

Ho^ <0ik)dLJ'^,. + j E dLJ^p.. (2) 

The hopping matrix K in the t2g basis is given by 

/ e(t,l,2,3) e(i',3,l,2) e(t',2,3,l)\ 
K'^(A;)= £(t',3,l,2) e(i,2,3,l) e(i',l,2,3) , 
Ve(i',2,3,l) e(i',l,2,3) £(t,3,l,2) / 

(3) 

with (1,2,3) = (fci,fc2,/c3), fci = V3kx/2~ky/2, fca = 
ky, kz = —ki — /c2, and e(t, 1,2,3) = 2iiC0sfci -t- 
2t2(cosA:2-l-cosA:3)-|-2i3 cos (fc2 — fc3)-l-2<4[cos (fcs — ki) + 
cos (fci - ^2)] + 2t5 cos (2fci) + 2t6[cos (2A:2) + cos (2A:3)] + 
■ ■ ■. The (t,t') denote the (intra, inter)-orbital hopping. 

FigCl shows the fitting of the tight-binding dispersions 
obtained by diagonalizing Eq. (01 to the LDA band struc- 
ture at a:=l/3 We note that the fit with up to 
third-nearest-neighbor (NN) hopping or more describes 
the LDA bands quite well. On the other hand, the tight- 
binding model cannot reproduce completely the LDA dis- 
persions even with up to eighth-NN hopping. The dis- 
crepancy is most pronounced along the M-K direction 
where the two e'g bands cross in the tight-binding fit 
(Fig. la). Similar disagreement can be traced back to 
the previous tight-binding fits 0, 0, 0| . We believe 
the difficulty arises from the hopping path via the O 2s 
and 2p orbitals. Nevertheless, the tight-binding model 
works very well at low energies near the Fermi level. The 
FS consists of a cylindrical sheet around the F-point and 
six hole pockets near the K-points as shown in Fig. Ib-d. 
The central FS has a dominant aig character while the six 
FS pockets are mainly of the e'g character. The hopping 
integrals obtained from the fit with up to third-NN are 
t=(-44.6, -9.0, 36.2, 5.9, 57.9, 36.7)meV and i'=(-157.8, 
-30.2, 37.1, 9.2, -11.9, -21.0)meV. The crystal-field split- 
ting A is chosen to be O.OleV. In the rest of the paper. 




FIG. 1: Tight-binding fits to the LDA band structure at a; = 
1/3. (a) The fitted band dispersions with up to 3rd, 5th 
and 8th NN hopping. The corresponding Fermi surfaces are 
plotted in (b), (c), and (d) respectively. 

we use these parameters for Hq. Our results are insen- 
sitive to these values provided that they provide a good 
fit of the LDA band structure near the FS. Note that 
although holes are evenly distributed among the three 
t2g orbitals, in the aig and e'g basis (hereafter referred 
to as the {a} basis), the hole occupations are 0.123 {e'g) 
and 0.421 {aig) respectively. Despite its higher orbital 
energy, the aig hole orbital has a higher occupation due 
to its larger bandwidth. 

The correlation effects are described by the multi- 
orbital Hubbard model H = Hq + Hj, where Hq is the 
tight-binding Hamiltonian in Eq.|(21), Hj represents the 
local Coulomb repulsion U (intra-orbital) and U' (inter- 
orbital) and Hund's rule coupling J h ■ F or t2g orbitals. 
Hi has been shown to take the form |lj| 

Hi = U'Yn^a^n.^al + {W - ^Jh) ^ n.ianii3 (4) 

i,Q. i,a'>j3 

with U' = U — 2Jii- Here fiia and Sia are the density 
and the spin operators in the {a} basis where the tight- 
binding part Hq is 

k,(T af3 k.a,(T 

Here A^ = —A/3, 2A/3 for the e'g and aig orbitals re- 
spectively. The hopping matrix Ti"'{k) — 0-^K'^(fc)0, 
with O the orthogonal rotation from the t2g to the {a} 
basis. Hi is identical in these two bases. The hierarchy 
of the interaction strength is U > U' > Jh > 0- 

We first study the effects of interactions in the HF 
theory in the orbital sector. In the paramagnetic phase, 
the interacting Hamiltonian is given by, 

k,a k,a^l3 
+ {U~2Ks) Y ["a/3«L«fe/3--^]. 



r KM r 

FIG. 2: HF results for U = S.OeV at x = 1/3. The band 
dispersions are shown for ri — 1/3, 1/2, 2/3. 

where n„/3 = i'^/Ns)Y.k,aKf3ik), Kpi^) = (aL^afc/s-r), 
Tia = riaa, and C/^g = U' — Jh/2. In essence, this HF 
analysis, also discussed in Ref. |l2|, is equivalent to the 
LDA+U theory . Since we are interested in the orbital 
dependent corrections, we have not displayed in Eq. ® 
the double-counting term which corrects for the energy 
already included in the LDA, because it depends only on 
the total density. 

Note that the HF theory is basis independent. It is 
thus convenient to stay in the {a} basis where the local 
density matrix, and thus the HF self-energies are diag- 
onal in orbital space, i.e. na^p = 0. Then the average 
energy has contributions only from the first two terms 
in Eq. ijHJi. Recall that most of the holes reside in the 
aig orbitals. If L/^g < U/2, it is favorable to increase 
the population of the orbitals to take advantage of the 
smaller inter-orbital repulsion U^g. On the other hand 
for U^ff > U/2, the tendency is to empty the e'g orbitals 
in favor of aig occupation. The crucial factor of 1/2 in 
front of U comes from the fact that due to exchange, 
intra-orbital repulsion operates only between holes with 
opposite spin, whereas both spins contribute to C/^g. 

We proceed to calculate the HF self-energy in terms 
of the average n = (riaig + 2ne' )/3 = (l — a;)/3 and the 
difference S = (j^aig — Tie' )/3 between the hole occupation 
of the aig and e'g orbitals, 

= lnU{l + 4v)+SUiv~l/2) (7) 
= lnU{l + Arj)-2SU{v-l/2) (8) 

where rj = U^g/U is the relative strength of the inter- 
orbital interaction. The interaction effect in the param- 
agnetic HF theory is to simply shift the atomic levels by 
and S^j^ respectively, resulting in a renormalization 
of the atomic level spacing A' = —35U{ri — 1/2). As ex- 
pected, the direction of the charge transfer depends on 
the ratio rj. Since the majority of the holes resides in the 
aig orbital in the noninteracting limit, S > 0. Thus, for 
T] < 1/2, the level splitting renormalizes upward. A' > 0, 
and interactions induce a transfer of carriers from the 
aig to the e'g orbital. The self-consistent HF results are 
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FIG. 3: The doping x dependence of the level spacing and 
bandwidth renormalization for the aig and e'g orbitals. 

shown in Fig. 2 at x = 1/3 for U = 3eV which is close 
to the value (3.7 eV) estimated from the LDA |13]. For 
T] = 1/3, the size of the hole pockets indeed becomes 
larger than that of the noninteracting/LDA ones. At 
r] — 1/2, the self-energy corrections are equal among the 
orbitals and the noninteracting (LDA) band dispersions 
are unchanged as shown in Fig. 2. When rj > 1/2, i.e. 
for U'/U > 3/5 or Jr/U < 1/5 which is reasonable for 
the cobaltates , the level splitting renormalizes down- 
ward, A' < 0, triggering a transfer of holes from the e'g to 
the aig orbital. The six FS pockets continues to shrink 
as the e'g band sinks with increasing rj and disappears 
beyond a critical ratio r]c, as shown in Fig. 2 for 77 = 2/3. 
We find that r]c{U = 3.0eV) ~ 5/8. 

The HF analysis shows that the disappearance of the 
six FS pockets near the K-points is the physics of large 
U and U' compared with Jh- In this case, the HF the- 
ory itself becomes unreliable. Moreover, the localization 
tendency leading to the bandwidth reduction cannot be 
captured by the HF theory. It is therefore instructive 
to study the problem in the strong-coupling limit by 
projecting out the states of double-occupation prohib- 
ited by the large on-site Coulomb repulsions. This can 
be achieved by the Gutzwiller projection {^{{ria})) = 
^Gl^o({"-a})), where Pq is the projection operator that 
operates on the noninteracting state |\l/o({nQ})) in a 
given orbital occupation scheme {ric}. It removes the 
double-occupancy by electrons from both the same and 
different orbitals. This variational procedure is most con- 
veniently implemented in the Gutzwiller approximation 
where the effect of projection is taken into account by 
the statistical weighting factor multiplying the quantum 
coherent state . Specifically, we approximate the hop- 
ping term by 

{^al„a,f3.\^)^gf{^o\a\^,a,p,\^^), (9) 

where the Gutzwiller renormalization factor gt is given 
by the ratio of the probabilities in the hopping process 
in the projected j^I') and the unprojected |5'o)- We find. 
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It is important to note that in a multi-orbital system 
g'^^ depends on the occupation of the orbitals connected 
by the hopping integral as seen in the denominator in 
Eq. H1U|) . The latter originates from the Pauli princi- 
ple. It compensates for the effects of "Pauli-blocking" 
of double-occupation by electrons in the same quantum 
states which already operate in the free fermion term on 
the RHS of Eq. while the numerator x in Eq. 1(1(1}) 
describes the "Coulomb-blocking" due to the large on- 
site U and U' . It turns out that the denominator is 
crucial for carrier transfer and orbital-polarization in 
the strong-coupling limit. In the uniform paramagnetic 
phase, g^^ = 2x / {2 — na){2 — np) . The orbital occu- 
pations are variational parameters determined by mini- 
mizing the ground state energy of the Hamiltonian, 

HgW = XI 9t'^'^al3{k)al^aakl3a + XI ^o.al^aO'iac, 

+ ^eiaC^aj^^aiaa - n^a), (11) 

where Sia are the Lagrange multipliers enforcing the oc- 
cupation n.ia — J^ai^icra'^i'^a) ■ They are determined by 
the self-consistency equation 

= = ^^i- E 5f KS^(aL,afe^.). (12) 

The RHS of this equation is the derivative of the kinetic 
energy of band a with respect to tIq,, and can be under- 
stood by the following argument. The transfer of a hole 
from band a to /3 causes the respective bandwidths to de- 
crease and increase by 0{1/Ns). However, the kinetic en- 
ergies of the occupied states in each band are changed by 
order unity, and this energy difference must be reflected 
in the equilibrium condition by an energy shift Eq.. Thus 
in contrast to HF theory, there is no energy cost pro- 
portional to U. Instead, both the band-narrowing and 
the renormalization of the level spacing A' — Eaig — Ee' 
contribute to the redistribution of holes among the or- 
bitals. In Fig. 3, we show the self-consistently determined 
band-narrowing factor Zq, = g^" and the renormalized 
level spacing. For orbitals with a larger hole occupation, 
the bandwidth reduction is smaller and the renormal- 
ized band energy is lower, resulting in the transfer of 
more holes into these bands. The combined effects cause 
the holes to move out of the e'g band into the aig band. 
The calculated band dispersions and the FS topology at 
X = 0.3, 0.5, 0.7 are shown in Fig^ in the strong cou- 
pling theory. The corresponding non-interacting case is 
also shown for comparison. The six hole pockets are com- 
pletely absent at all levels of sodium doping due to strong 
correlation, leaving a single hexagonal Fermi surface cen- 
tered around the F point satisfying the Luttinger theo- 
rem. This, as well as the band-narrowing due to strong 
Coulomb repulsion, is in very good agreement with the 
photoemission experiments |j,|5(. 
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FIG. 4: The band dispersions and the Fermi surfaces in the 
strong coupling limit (red solid lines) for doping x — 0.3, 0.5 
and 0.7. The noninteracting dispersions (blue dashed lines) 
are also plotted for comparison. 

In conclusion, we have shown that strong correlation 
plays an important role in understanding the electronic 
structure of Naa;Co02. It pushes the e'g band below the 
Fermi level, leading to an orbital polarized state with a 
single hole-like FS. The absence of the small FS pockets, 
which would have contributed significantly to the density 
of states in band structure calculations, further suggests 
that the large mass enhancement observed in the specific 
heat measurement ^3 due to strong correlation. 
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